Reliability test for the experimental results of electric quadrupole hyperfine structure 
constants and new assessment of nuclear quadrupole moments in i35,i37g^ 
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We investigate electric quadrupole hyperfine structure constant (B) results in ^^^Ba"*" and 
^■^^Ba^ and compare their ratios with the reported measurements of the ratio between the nuclear 
quadrupole moment (Q) values of these isotopes. We carry out confidence tests for the reported 
experimental B values and calculations oi B /Q from the present work. Inconsistencies in the exper- 
imental B values are observed in both the isotopes from different experiments performed using the 
same techniques. The present calculations are carried out using an all order relativistic many-body 
theory considering only single and double excitations in the coupled-cluster ansatz. After a detailed 
analysis of the results, the values of Q we obtain for ^^^Ba and ^^'^Ba are 0.153(2)& and 0.236(3)&, 
respectively, which differ by about 4% from the currently referred precise values. 
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The effective ellipsoidal nuclear charge distribution of 
an isotope in the nuclear ground state is generally quan- 
tified by the nuclear electric quadrupole moment (Q). 
The nuclear shell model assumes that the charge distribu- 
tion is spherically symmetric for a doubly-magic closed- 
shell nucleus, and therefore its quadrupole moment van- 
ishes; otherwise it will have an intrinsic non-zero Q value. 
Thus, accurate evaluations of Q for different isotopes 
could probe the shell theory (e.g. see [ij) and therefore 
the investigation of this quantity has been of consider- 
able interest to nuclear, _atomic and molecular physicists 
for the last six decades 



Although the value of Q for nuclear ground states can- 
not be measured directly, it is possible to measure their 
ratios for different isotopes [a-fllj. Also, for an atomic 
state of a particular isotope, Q can be determined by 
combining the measured electric quadrupole hyperfine 
structure constant [B) with a calculation oi B/Q for that 
isotope. To obtain an accurate value of Q by this proce- 
dure, both the measurement and the calculation have to 
be performed to high precision. Experimental results are 
generally considered to be more reliable than the calcu- 
lated results for many-electron atoms owing to the fact 
that ab initio evaluation of various physical quantities 
using many-body methods involve a number of approxi- 
mations at different stages of the calculations. However, 
measurements from different experiments are not always 
in agreement and, in certain cases, large discrepancies 
between measurements and accurate calculations have 
been noted [H, [ll] . It is therefore essential to scrutinize 
the accuracies of both the experimental B and calculated 
B/Q results when determining values of Q from them. 

Studies of hyperfine structure constants in singly ion- 



ized barium (Ba^) can provide valuable information 
about the wave functions in the nuclear region and the 
role of the electron correlation effects whose knowledge 
are relevant for the proposed parity nonconservation ex- 
periments [l^ . It is also important for the extraction of 
higher order nuclear moments [Tsj . Furthermore, studies 
of the hyperfine structure are also interesting for the in- 
vestigation of the variation of the fine structure constant 
and for astrophysical studies JL7J . 

In this Rapid Communication, we investigate the con- 
sistency between measured values of B and compare 
these to calculations oi B/Q for the 135 and 137 iso- 
topes of Ba+. This provides us with new estimates of 
the quadrupole moments for these two isotopes that dif- 
fer from previous estimates by approximately 4%. In 
addition we show that there are significant discrepancies 
in the experimentally measured values of B. 

The nuclear electric quadrupole hyperfine structure 
constant B for an atomic state with angular momentum 
J > 1/2 is defined as fis!] 

= QR'iJ) (1) 

The operator T^^^ = t^^^ 

is a second rank tensor whose 
reduced matrix element in terms of the single particle 
orbitals is given by 
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where Ki and ji are the angular momentum quantum 
numbers of i*^^ Dirac orbital with large Pi and small Qi 
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TABLE I: Calculated results of = B/Q in MHz/6 and A factors in ^^^Ba+ and "^Ba+ from the present work. Estimated 
uncertainties from the calculations are given in the parentheses. 



Atomic 


"^Ba 


^■^'■Ba 


State 


^137 

^137 


[5p^] 6p3/2 


[5p''] 5d3/2 


[5p''] 5d5/2 


A 137 
^135 


-Rl35 


[5pS] 5d3/2 
[5/] 5d5/2 


379.52(4.0) 
188.69(2.5) 
253.76(3.5) 


0.0 
1.01(5) 
0.50(4) 


-0.50(1) 

0.0 
-0.26(2) 


-0.33(2) 
0.35(3) 
0.0 


-1.2 X lO"'^ 
3.0 X 10"'^ 
-4.0 X lO"'^ 


379.52(4.0) 
188.69(2.5) 
253.76(3.5) 



radial components. In Eq. ([T]), we have defined the ratio 
R'^{J) = B/Q which is a function of J in the electronic 
coordinate and can be determined directly from theoret- 
ical calculations. 

Within first order perturbation theory, the hyperfine 
structure constants for a particular state are linearly re- 
lated to measurements of the hyperfine splittings (see, 
for example, [l9j and references therein). Thus, provided 
higher order terms can be neglected, the hyperfine struc- 
ture constant, B, can be determined directly from exper- 
iments. Measured values of B and calculated values of 
can then be used to determine Q. Since the ratio of 
Q values between isotopes can be measured directly, a 
useful consistency check on both measurements and cal- 
culations is to consider ratios between different isotopes 
or different states within the same isotope. In general we 
have the relationship 
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where we have introduced — {R\ — i?2)/^2 ^^"^ ^^'^ 
subscripts 1 and 2 denote the relevant parameter for dif- 
ferent states and/or isotopes. The superscripts expt and 
theor are used to emphasize the fact that the correspond- 
ing term can be obtained from experimental measure- 
ments and theoretical calculations, respectively. The pa- 
rameter A2 is analogous to the hyperfine anomaly that 
has been studied extensively for the magnetic dipole hy- 
perfine structure constant {A) in a number of atomic sys- 
tems m, [2I], but has not been studied for B. Indeed, 
the above expression can be generalized to all hyperfine 
structure constants. 

In Table HI we present calculated values of and A 
among the 6^3/2, 6^3/2 and 6^5/2 states for ^^^Ba+ and 
^^^Ba+. Atomic states that have been considered here 
have a common closed core [5p^] and one electron in the 
valence orbital, denoted by k, in the different quantum 
states in Ba+. The wave functions of these states with 
the Dirac-Coulomb (DC) Hamiltonian are calculated us- 
ing the relativistic coupled-cluster (RCC) method by ex- 
pressing them in a form (2^ 



where |<I>fc) is the Dirac-Fock (DF) reference state corre- 
sponding to the valence electron k. T and Sk are the op- 
erators that include correlations among the core electrons 
and the valence and the core electrons respectively. We 
consider all possible single and double excitations (known 
as CCSD method) along with the leading triple excita- 
tions (known as CCSD(T) method) involving the valence 
orbital. 

The finite charge distribution of the nucleus in the 
different isotopes is modeled by taking a two-parameter 
Fermi-nuclear charge distribution 



Po 



(5) 



where po is the density for the point nuclei, c and a 
are the half-charge radius and skin thickness of the nu- 
cleus. These parameters are chosen as a = 2.3/4(Zn3) 

and c = 



5 2 
3 rms 



la^TT^, where 



(4) 



r rms is the root mean 
square radius of the atomic nucleus that is given in the 
Fermi model as r^ms = 0.836^^/'^ -I- 0.57 in fermi (/m) 
with atomic mass A. To reduce the uncertainty in our 
calculation for different charge distributions, we consider 
the optimal values of rrms as 4.8273 and 4.8326 frn for 
^■^'^Ba+ and ^"^^Ba"'", respectively, as tabulated in [25]. 

The uncertainties in the calculations presented in Ta- 
ble U are estimated from the finite size of bases, approxi- 
mations in the level of calculations and higher order rela- 
tivistic corrections. As seen in the table, the Ajgy factors 
with respect to different states are large and the associ- 
ated uncertainties are also significant. However, the Ajg^ 
factor for a state is below 10~^ and it can only play a role 
when the accuracies of the quantities of interest are com- 
parable to these values. This is due to the fact that the 
wave function for a particular state does not depend sig- 
nificantly on the isotope and we can expect this to be true 
for the neutral atom also. Thus, from Eq. ([3]), we see 
that the ratio -B137/-B135 gives the ratio Qi37/(3i35 which 
is independent of the state of interest. This provides a 
direct consistency check of the experimentally measured 
ratios of B and Q. 

In Table |TT] we present reported measurements of the 
ratio Q137/Q135. The last four are measured directly us- 
ing nuclear magnetic resonance techniques whereas the 
first four are extracted from hyperfine measurements. 
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TABLE II: Reported ratios of Q values from various experi- 
ments. Q values are given in h. 



Method 



1.64(16) 
1.531(26) 



0.245(4)/0.160(3) = 1.531(38) 
0.228(24)/0.146(16) 
0.246(2)/0.150(15) : 
0.248(3) /0.162(2) = 
1.543(3) 
1.537(2) 
1.49(10) 
1.538485(95) 
1.5426(62) 



Fast beam spectroscopy of 
6p3/2 state in Ba^ [26] 
1.56(24) Level-crossing spectroscopy 
of 5d6p ^Pi state in Ba [27] 
Fast beam spectroscopy of 
5ci3/2 state in Ba"*" [28] 
Non-relativistic calculation 
using CCSD method [29| 
Proton magnetic resonance 
resonance in BaBr2 [8] 
Level-crossing spectroscopy 
of 5d6p ^Pi state of Ba [3ff] 
NMR spectroscopy of 
BaCl2 [31] 

NMR and NQR spectroscopy 
ofBaCl2 [32] 
Microwave spectroscopy 
of BaO [33] 



TABLE III: Brief summary of the measurements of B (in 
MHz) in different states of ^"^^Ba and ^^^Ba and in their singly 
charged ions from the measurements carried out by same ex- 
periments. Ratios of these quantities are also given to com- 
pare them with the ratio of Q values listed in Table [III 





State 




oa 


135 -R., 

oa 
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Ratio 
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47 ^Q(^(^f,'\'' 

^ 1 .Oi/Ul-LUJ 


30 Sni C941'' 


1 53Sfi('13"> 


[5p®]6s5d 
[5p®]6s5d 


^i?2 


59.564(14)" 
17.890(3)'' 


38.710(15)" 
11.642(4)'' 


1.5387(68) 




1.53668(59) 


[5/]6s6p 


'Pi 


50.09(21)= 
49.7(5)'* 


34.01(22)"= 


1.473(11) 






32.5(4)"* 


1.529(24) 


[5p«]6s6p 


'Pi 


-41.61(2)"= 


-27.08(2)"= 


1.5366(14) 


[5/]6p5d 


'Di 


-3.5(6)^ 


-2.4(10)-^ 


1.46(66) 


[5/]6p5d 




45.3(10)-'' 


29.0(10)-'^ 


1.562(64) 


[5p>p5d 


'F, 


83.0(14)-'' 


53.9(14/ 


1.540(48) 


[5/]6p5d 


-'F, 


111.0(13)^ 


71.1(13/ 


1.561(34) 


[5/]6s7p 


'Pi 


16.6(6)3 


8.8(6)3 


1.89(15) 



Ba II 




5 10 15 20 25 30 35 

Results given sequentially from Tables II and III 



FIG. 1: Plot showing {x — xq) /Sx with x representing ratio of 
Q values for ^^^Ba and ^^^Ba obtained from different works 
given sequentially in Tables HH and HiH (taken in X-axis), xo 
denoting the most accurately reported ratio 1.538485 in Ref. 
[3^ 1 and Sx correspond to errors associated with the respective 
X values. 



In Table lllli we present the reported measurements of 
B in both atoms and singly charged ions of ^^^Ba and 
^■^^Ba, along with the ratio Bi^j / The uncertainties 
given are assumed to represent one standard deviation. 
To highlight the variation in the results, in Fig. [TJ we 
plot the deviation of each ratio relative to its respective 
uncertainty from the most accurately reported value of 
1.538485. From this plot we see there are four experimen- 
tal values in which the deviation is more than 5 standard 
deviations. Thus these reported values should be treated 
as questionable. We note that a minimization of the 



[5p']6p3/2 


92.5(2)'' 


59.0(1)" 


1.5678(43) 




95.0(37)* 


59.4(23)* 


1.599(88) 




89.7(15)^' 
91.6(12)*= 


55.3(34)^ 
62.1(12/ 


1.62(10) 




1.475(34) 




92.5(4)' 


60.4(4)' 


1.531(12) 




97.1(93)' 


63(9)' 


1.54(27) 


[5p'^]5d3/2 


47.5(13)^' 


33.2(24)^ 


1.43(11) 




44.5408(17)'" 


28.9528(25)'" 


1.53839(15) 


[5p*]5d5/2 


80.7(10)^ 


45.2(10)^ 


1.785(45) 




59.533(43)*" 


38.692(10)™ 


1.5386(12) 


References: 


""[34J; "[35]; ^[36]; 
"[40]; '[4^; ^[42]; 


'[37]; "=[30]; 
'[43]; 


^[38j; ^[39]; 




'"[28] 



remaining ratios gives a reduced of 1.02 and an esti- 
mated ratio of 1.538425(78) consistent with the reported 
value of 1.538485(95). We therefore conclude that the 
reported value of 1.538485(95) is a reliable estimate of 
the ratio Qisi/Qub- 

Having established a reliable value of Q137/Q135, we 
now proceed to assess the validity of the calculations pre- 
sented in Table [D To do this we consider experimentally 
measured B values for different states within the same 
isotope, specifically ^'^^Ba+. From Eq. ([3]), it is seen 
that the ratio of these B values give a direct measure 
of the calculated value of A\lj. Thus comparison of the 
measured ratio of B values with 1 -I- A^j provides a con- 
sistency check for our calculated i?"^ values. Taking the 
most accurately reported values of B in Table IIIl namely 
those for the ^3/2 and the 1^5/2 states from 28|, we de- 
duce a value of 1 -I- AJ3Y — 1.3366(1) consistent with our 
calculated value of 1.35(3). Thus we may now reliably 
use the measured values of B with our calculated values 
of i?"= to deduce Q values for the two isotopes. 
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TABLE IV: New B results for "■^Ba+ and "5Ba+ (in MHz) 
obtained combining the new Q values with our calculations 
of W factors . 



State 


i^^Ba 


""^Ba 


[5p'']6p3/2 


89.57(1.48) 


58.07(1.0) 


[5p'']5rf3/2 


44.53(82) 


28.87(54) 


[5/]5d5/2 


59.89(1.0) 


38.83(74) 



To obtain the greatest accuracy, we determine the Q 
value using data from the ^3/2 state of ^^''Ba+. Recently 
the B value for this state has been measured to very high 
accuracy 44.5387936(10) MHz fT^I consistent with the 
value given in Table Hill Combined with our calculated 
value given in Table U we obtain Q for ^^^Ba of 0.236(3) 
h. Further, from the ratio of Q values, 1.538485(95), 
we then deduce a value of Q for ^^^Ba of 0.153(2) 6. 
Comparing these values with other reported values of Q 
in Table HIl we find significant differences in the results. 
Recently, the value of Q for ^^^Ba was evaluated by us to 
be 0.246(1) b This value was determined using the 
experimental result for B of the [5p^] 6^3/2 state from [40|. 
However, as we have shown here, those measurements 
have significantly differ from the other measurements. 

For completeness, we can use the newly obtained Q val- 
ues and our calculations of the i?^ factors given in Table 
m to determine the B values of the [5p^] 6^3/2, [5p^] 5^3/2 
and [5p^] 5^5/2 states in ^^^Ba+ and ^^^Ba+. The results 
are given in Table IIVI and can be compared with the re- 



sults given in Table IHII As seen in Table IIVI the present 
estimated values are certainly not the most accurate val- 
ues due to the large uncertainties in our calculated 
factors, but they give consistent results for both the iso- 
topes and states considered. For more precise values of B 
and to determine more precise values of Q it is necessary 
to carry out more accurate theoretical calculations. 

In summary, we have carried out an extensive consis- 
tency check of reported experimental values of the hyper- 
fine structure constants over a range of different states 
for the 135 and 137 isotopes of Barium. Using the rela- 
tivistic coupled-cluster method, we have carried out the 
corresponding calculations and investigated ratios of the 
quadrupole moments between ^■^''Ba+ and ^'^^Ba'''. From 
our analysis, wc have found some inconsistencies in the 
reported measurements. The reason for such inconsis- 
tencies is not clear to us; one possibility is the neglect of 
higher order nuclear moments or higher order correction 
terms to the hyperfine interaction. We have also esti- 
mated Q values for both the ^"^^Ba^ and ^■^^Ba^ isotopes 
from results that were shown to be consistent. We pro- 
pose more experimental and theoretical studies of these 
quantities to ascertain our results. Nonetheless, our pro- 
posed analysis for the consistency check of the experi- 
mental results will also be useful in the future studies of 
the above properties. 
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tory, Ahmedabad. 



[1] H. Haas, Ch. Stenzel, H. Grawe, H. E. Mahnke and K. 

H. Maier, Nuclei Off the Line of Stability, ed. by R. A. 

Meyer and D. S. Brenne, ACS Symposium Series 324, 

Chap. 60, pp. 399, 1986. 
[2] N. F. Ramsey, Ann. Rev. Nuc. Sc. 1, 97 (1952). 
[3] B. T. Feld, Ann. Rev. Nuc. Sc. 2, 239 (1953). 
[4] P. Dufek, P. Blaha and K. Schwarz, Phys. Rev. Lett. 75, 

3545 (1995). 

[5] Y. Tanaka, R. M. Steffen, E. B. Shera, W. Renter, M. 
V. Hoehn and J. D. Zumbro, Phys. Rev. Lett. 51, 1633 
(1983). 

[6] T. Minamisono, T. Ohtsubo, I. Minami, S. Fukuda, A. 
Kitagawa, M. Fukuda, K. Matsuta, Y. Nojiri, S. Takeda, 
H. Sagawa and H. Kitagawa, Phys. Rev. Lett. 69, 2058 
(1992). 

[7] G. Neyens, Rep. Prog. Phys. 66, 633 (2003). 
[8] R. M. Wilhams and J. W. McCrath, Phys. Rev. 127, 184 
(1962). 

[9] R. M. Sternheimer and R. F. Peierls, Phys. Rev. A 3, 
837 (1971). 

[10] A. Voss et al, J. Phys. G 38, 075102 (2011). 

[11] M. Keim, U. Georg, A. Klein, R. Neugart, M. Neuroth, S. 
Wilbert, P. Lievens, L. Vermeeren, B.A. Brown, and the 
ISOLDE CoUaboration, Euro. Phys. J. A 8, 31 (2000). 

[12] B. K. Sahoo and B. P. Das, Phys. Rev. A 86, 022506 



(2012). 

[13] B. K. Sahoo, Md. R. Islam, B. P. Das, R. K. Chaudhuri 
and D. Mukherjee, Phys. Rev. A 74, 062504 (2006). 

[14] B. K. Sahoo, P. Mandal and M. Mukherjee, Phys. Rev. 
A 83, 030502(R) (2011) and references therein. 

[15] N. C. Lewty, B. L. Chuah, R. Cazan, B. K. Sahoo and 
M. D. Barrett, Optics Express 20, 21379 (2012). 

[16] S. Karshenboim, Can. J. Phys. 78, 639 (2000). 

[17] L. Mashonkina, T. Gehren and I. Bikmaev, Astro. As- 
trophys. 343, 519 (1999). 

[18] C. Schwartz, Phys. Rev. 97, 380 (1955). 

[19] E. Arimondo, M. Inguscio and P. Violino, Rev. Mod. 
Phys. 49, 31 (1977). 

[20] N. J. Stone, Table of Nuclear Magnetic Dipole and Elec- 
tric Quadruopole Moments, IAEA Nuclear Data Section, 
Vienna International Centre, 1400 Vienna, Austria, April 
2011. 

[21] P. Pyykko, Mol. Physics 106, 1965 (2008). 

[22] S. Buttgenbach, Hyperfine Int. 20, 1 (1984). 

[23] J. R. Persson, arXiv:110:5991 (unpublished). 

[24] B. K. Sahoo, G. Gopakumar, R. K. Chaudhuri, B. P. 

Das, H. Merlitz, U. S. Mahapatra and D. Mukherjee, 

Phys. Rev. A 68, 040501 (R) (2003). 
[25] I. Angeh, At. Data and Nuc. Data Tables 87, 185 (2004). 
[26] K. Wendt, S. A. Ahmad, C. Ekstrom, W. Klempt, R. 



5 



Neugart and E. W. Otten, Z. Phys. A 329, 407 (1988). 
[27] I. -J. Ma and G. zu Putlitz, Z. Physik A 277, 107 (1976). 
[28] R. E. Silvorans, G. Borghs, P. Dc Bisschop and M. Van 

Hove, Phys. Rev. A 33, 2117 (1986). 
[29] A.-M. Martensson-Pendrill and A. Ynnerman, J. Phys. 

B 25, L551 (1992). 
[30] G. zu Puthtz, Ann. Physik (Leipzig) 11, 248 (1963). 
[31] H. Kriiger, O. Lutz and H. Oehler, Phys. Letts. 62A, 131 

(1977). 

[32] O. Lutz and H. Oehler, Z. Physik A 288, 11 (1978). 
[33] C. E. Blom, H. G. Hedderich, F. J. Lovas, R. D. Suenram 

and A. G. Maki, J. Mol. Spectroscopy 152, 109 (1992). 
[34] S. G. SchmelHng, Phys. Rev. A 9, 1097 (1974). 
[35] M. Gustavsson, G. Olsson and A. Rosen, Z. Phys. A 290, 

231 (1979). 

[36] W. A. van Wijngaarden and J. Li, Can. J. Phys. 73, 484 
(1995). 



[37] G. Nowicki, K. Bckk, S. Goring, A. Hanscr, H. Rebel and 

G. Schatz, Phys. Rev. C 18, 2369 (1978). 
[38] P. Grundevik, M. Gustavsson, G. Olsson and T. Olsson, 

Z. Phys. 312, 1 (1983). 
[39] E. R. Eliel, W. Hogervorst, T. Olsson and L. R. Pendrill, 

Z. Phys. A 311, 1 (1983). 
[40] P. Villemoes, A. Arnesen, F. Heijkcnskjold and A. Wbn- 

strom, J. Phys. B 26, 4289 (1993). 
[41] K. Wcndt, S. A. Ahmad, F. Buchinaer, A. C. Mueller, R. 

Neugart and E. W. Otten, Z. Phvs. A 318, 125 (1984). 
[42] R. E. Silverans, G. Borghs, G, Dumont and J. -M. Van 

den Cruyce, Z. Phys. 295, 311 (1980). 
[43] H. Winter and M. Gilliard, J. Phys. B 10, 2739 (1977). 
[44] H. J. Andra, Beam-Foil Spectroscopy, ed. by I. A. Sellin 

and D. J. Pegg (New York: Plenum) p. 835 (1976). 
[45] B. K. Sahoo, Phys. Rev. A 74, 020501(R) (2006). 



